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Abstract — Via collaborative beamforming, nodes in a wireless 
network are able to transmit a common message over long 
distances in an energy efficient fashion. However, the process of 
making available the same message to all collaborating nodes in- 
troduces delays. In this paper, a MAC-PHY cross-layer scheme is 
proposed that enables collaborative beamforming at significantly 
reduced collaboration overhead. It consists of two phases. In the 
first phase, nodes transmit locally in a random access time-slotted 
fashion. Simultaneous transmissions from multiple source nodes 
are viewed as linear mixtures of all transmitted packets. In the 
second phase, a set of collaborating nodes, acting as a distributed 
antenna system, beamform the received analog waveform to one 
or more faraway destinations. This step requires multiplication 
of the received analog waveform by a complex weight, which is 
independently computed by each cooperating node, and which 
allows packets bound to the same destination to add coherently 
at the destination node. Assuming that each node has access 
to location information, the proposed scheme can achieve high 
throughput, which in certain cases exceeds one. An analysis of the 
symbol error probability corresponding to the proposed scheme 
is provided. 

I. Introduction 

Transmission over long distances often requires significant 
amounts of energy in order to overcome attenuation. Energy 
is usually a scarce commodity in wireless ad hoc networks, 
as nodes typically operate on batteries, which in many cases 
are difficult to replace or recharge. Thus, energy-efficient 
schemes for long-distance transmission in wireless networks 
have recently been of much interest. In some such situations, 
multihop may be a preferred solution. However, there are sev- 
eral challenges in transmitting real-time services over multiple 
hops. For example, the traditional CSMA/CA based medium 
access control (MAC) for avoiding collisions does not work 
well in a multihop scenario because transmitters are often 
out of reach of other nodes' sensing ranges. Thus, packets 
traveling across the network experience interference and a 
large number of collisions, which introduce long delays. Also, 
multihop networks require a high node density which makes 
routing difficult and affects the reliability of links (TJ. 

Recently, a collaborative beamforming technique was pro- 
posed in [3 1, in which randomly distributed nodes in a network 
cluster form an antenna array and beamform data to a faraway 
destination without each node exceeding its power constraint. 
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Research under Grant N00014-07-1-0500. 



The destination receives data with high signal power. Beam- 
forming with antenna arrays is a well studied technology; it 
provides space-division multiple access (SDMA) which en- 
ables significant increases in communication rate. A challenge 
with implementing beamforming in ad hoc networks is that the 
geometry of the network may change dynamically. In 0, it 
was shown that randomly distributed nodes can achieve a nice 
average beampattern with a narrow main lobe and low side 
lobes. The directivity of the pattern increases as the number of 
collaborating nodes increases. Such an approach, when applied 
in the context of a multihop network reduces the number of 
hops needed, thereby reducing packet delays and improving 
throughput. However, to study network performance, one 
must take into account the information-sharing time that is 
required for node collaboration. If a time-division multiple- 
access (TDMA) scheme were to be employed, the information- 
sharing time would increase proportionally to the number of 
source nodes (i.e., the nodes having packets to transmit). 

In this paper we propose a scheme that is based on the idea 
of collaborative beamforming, and reduces the time required 
for information sharing. A preliminary version of the proposed 
scheme appeared in [4|. The work in this paper contains error 
analysis that provides insight into the performance of the pro- 
posed approach. The main idea is as follows. Different source 
nodes in the network are allowed to transmit simultaneously. 
Collaborating nodes receive linear mixtures of the transmitted 
packets. Subsequently, each collaborating node transmits a 
weighted version of its received signal. The weights are such 
that one or multiple beams are formed, each focusing on one 
destination node, and reinforcing the signal intended for a 
particular destination as compared to the other signals. Each 
collaborating node computes its weight based on the estimated 
channel coefficients between sources and itself. This scheme 
achieves higher throughput and lower delay with the cost of 
lower SINR as compared to [3|. In the preliminary version of 
this work [4], the analysis of interference at the receiving node 
was done asymptotically, i.e., as the number of collaborating 
nodes tends to infinity. Here we provide analytical expressions 
for symbol error probability (SEP) that directly depend on 
the number of collaborating nodes. The analysis shows how 
SEP is affected by transmission power, signal-to-noise ratio, 
number of simultaneously transmitting nodes and number of 
collaborating nodes. 



II. Background on collaborative beamforming 

For simplicity, let us assume that sources and destinations 
are coplanar. We index source nodes using a subscript i, with tj 
denoting the i-th node. At slot n, one source node t m needs 
to transmit the signal s m (n) to a faraway destination node 
q m . Suppose that set of N nodes, designated as collaborating 
nodes ci, . . . , cm, have access to s m (n). The locations of these 
collaborating nodes follow a uniform distribution over a disk 
of radius R. We denote the location of Cj in polar coordinates 
with respect to the origin of the disk by (rj, ipi). Let di m ((p m ), 
or simply di m , represent the distance between Cj and the 
destination q m , where <fi m is the azimuthal angle of q m with 
respect to the origin of the disk. d 0m (4> m ) or do m denotes the 
distance between the origin of the disk and q m , so the polar 
coordinates of q m are (do m , <fi m ). Moreover, let di{(j>) denote 
the distance between Cj and some receiving point with polar 
coordinate (do m ,0). The initial phases at the collaborating 
nodes are set to 

*i(^m) = -yC(^™), i = l,...,N . (1) 

This requires knowledge of distances (relative to wavelength 
A) between nodes and destination, and applies to the closed- 
loop case [3|. Alternatively, the initial phase of node i can 
be 

*j(<^m) = — n CO$((f> m - Ipi) (2) 

which requires knowledge of the node's position relative to 
some common reference point, and corresponds to the open- 
loop case [3|. In both cases synchronization is needed, which 
can be achieved via the use of the Global Positioning System 
(GPS). 

The path losses between collaborating nodes and destination 
are assumed to be identical for all nodes. The corresponding 
array factor given the collaborating nodes at radial coordinates 
r = [r%, Vn] and azimuthal coordinates ip = [-01, i/jn] 
at location with polar coordinate (do m > </>) is 

1 N 

F(0;m|r,-0) = ± ^ e i*«(*») e ** d «(*) . (3) 
i=l 

Under far-field assumptions, the array factor becomes J3] 

1 N 

f{4>m^) = -xy a <**»>« (4) 

1=1 

where a((j);<j> m ) = Att{R/\) sin(|(^ m - <j>)), and ^ = 
(rt/R) sm(ipi — \ {4> m + 4>))- The random variable z\ has the 
following probability density function (pdf): 

f Zi (z) = -y/l-z*, ~l<z<l. (5) 

7T 

Finally, the average beampattern can be expressed as |f3] 

Pav(0) = E z {\F(0\ Z )\ 2 } 

= i+fl-i) 2^1%H 2 (6) 
N \ NJ a(0;0 m ) 



where Ji(.) is the first-order Bessel function of the first kind. 
When plotted as a function of 4>, P av {4>) exhibits a main lobe 
around <j) m , and side lobes away from <p m . It equals one in 
the target direction, and the sidelobe level approaches 1/N as 
the angle moves away from the target direction. The statistical 
properties of the beampattern were analyzed in |3|, where it 
was shown that under ideal channel and system assumptions, 
directivity of order N can be achieved asymptotically with N 
sparsely distributed nodes. 

As we have noted, all of the collaborating nodes must 
have the same information to implement beamforming. Thus, 
the source nodes need to share their information symbols 
with all collaborating nodes in advance. If a TDM A scheme 
were to be employed, the information-sharing time would 
increase proportionally to the number of source nodes. In 
the following, we propose a novel scheme to reduce the 
information-sharing time and also allow nodes in the network 
to transmit simultaneously. 

III. The proposed scheme 

Here we refine the model of [3|, focusing more directly on 
the physical model for the signal, fading channel and noise. 
In addition to the above assumptions, we will further assume 
the following: 

1) The network is divided into clusters, so that nodes in 
a cluster can hear each other. In each cluster there is a 
node designated as the cluster-head (CH). Nodes in a 
cluster do not need to transmit their packets through the 
CH. 

2) A slotted packet system is considered, in which each 
packet requires one slot for its transmission. Perfect 
synchronization is assumed between nodes in the same 
cluster. 

3) Nodes transmit packets consisting of phase-shift keying 
(PSK) symbols each having the same power a 2 . Also, 
nodes operate under half-duplex mode, i.e., they cannot 
receive while they are transmitting. 

4) Communication takes place over flat fading channels. 
The channel gain during slot n between source i, and 
collaborating node Cj is denoted by ay(n). It does 
not change within one slot, but can change between 
slots. The channel gains are independent and identically 
distributed (i.i.d.) complex Gaussian random variables 
with zero means and variances a 2 across both time and 
space, i.e., a,ij(n) ~ CA/"(0, a 2 ). 

5) The complex baseband-equivalent channel gain between 
nodes q and q m is bi m e^ dim (6), where bi m is the 
path loss. The distances between collaborating nodes 
and destinations are much greater than the maximum 
distance between source and collaborating nodes. Thus, 
bi m is assumed to be identical for all collaborating nodes 
and equals the path loss between the origin of the disk 
and the destination, denoted by b m . 

Suppose that cluster C contains J nodes. During slot n, 
source nodes t±,...,tK need to communicate with nodes 
qi, . . . ,qK that belong to clusters C\, . . . , Ck, respectively. 



The azimuthal angle of destination is denoted by fa. 
The packet transmitted by node tj consists of L symbols 
Sj(n) = [sj(n;0), . . . , Sj(n; L — 1)]. Due to the broadcast 
nature of the wireless channel, non-source nodes in cluster C 
hear a collision, i.e., a linear combination of the transmitted 
symbols. More specifically, node Cj hears the signal 

K 

x i( n ) = /, a >ji( n ) s j( n ) +Wj(n) (7) 
i=i 

where w,(n) = [wi(n; 0), . . . , Wi(n; L — 1)] represents noise 
at the receiving node c, . The noise is assumed to be of zero 
mean and with co variance matrix ct^Il, where II is an L x L 
identity matrix. 

Once the CH establishes that there has been a transmission, 
it initiates a collaborative transmission period (CTP), by send- 
ing a control bit to all nodes, e.g., 1, via an error-free control 
channel. The CH will continue sending a 1 in the beginning 
of each subsequent slot until the CTP has been completed. 
The cluster nodes cannot transmit new packets until the CTP 
is over. 

Let q m denote the destination of s m (n). In slot n + m, to — 
1, . . . , K, each collaborating node Cj transmits the signal 

x;(n + to) = ^(n)fi m a* ml (n)e^^ (8) 

where \i m is a scalar used to adjust the transmit power and 
is the same for all collaborating nodes. fi m is of the order of 
1/N. 

Collaborating nodes need to know which are source nodes 
and then estimate the channel between all source nodes and 
themselves. One possible way to implement this is to use 
orthogonal IDs, as discussed in ED . 

Also, collaborating nodes require the knowledge of their 
initial phases. In closed-loop mode, each collaborating node 
can independently synchronize itself to a beacon sent from 
the destination and adjusts its initial phase to it [3|. In open- 
loop mode, each collaborating node needs to know its relative 
position from a predetermined reference point (e.g. the origin 
of the disk) within the cluster, which can be achieved by the 
use of GPS. To obtain initial phases, collaborating nodes also 
require knowledge of the azimuths of the destinations so that 
the beams can be steered toward desired directions, which may 
be broadcast by the CH via a control channel. 

Given the collaborating nodes at radial coordinates r = 
[n, ...,rjv] and azimuthal coordinates ip = [tpx, i/ipf], the 
received signal at an arbitrary location with polar coordinates 
(d Qm , fa), is 

N 

y(fam\r,ip) = b m ki(n + m)e i ^ + v(n + to) (9) 
i=i 

where v(n + to) represents noise at the receiver during slot 
n + to. The covariance matrix of v(n + to) equals (t^Il- 

It was shown in [4] that, as N — > oo and omitting the noise, 
y{fan\ rn\r, i/>) — * N n m b m cr 2 s m (n) . Thus, the destination 
node q m receives a scaled version of s m (n). The beamforming 



step is completed in K slots, reinforcing one source signal at 
a time. 

Assuming that all of the K source packets have distinct 
destinations at different resolvable directions, multiple beams 
can be formed in one slot, each beam focusing on one direction 
and reinforcing one source signal. In the rest of the paper, for 
simplicity we will consider only the case in which a single 
beam is formed during slot n + to, focusing on destination 
q m . The results obtained under this assumption can be readily 
extended to multiple simultaneous beams. 

Taking into account the assumptions on channels and noise, 
the average beampattern was derived in [4|. Defining the 
throughput, T, as the average number of packets that are 
successfully transmitted in a time slot, we showed in |4| that 
K/(l + K) < T < K/2, which could be greater than 1. 
Also, in Q, we showed that under a fixed transmit power, 
the average signal-to-interference plus noise ratio (SINR) is 
asymptotically (3' times less than that of 0, where f3' = 
K+1 + 

IV. Symbol Error Probability (SEP) 

In the following, for simplicity we omit the time index, and 
replace y, it;, Xi, s;, w; and v in the above equations by y, Xi, 
Xi, Si, Wi and v (i.e., with one of their samples) respectively. 

Our analysis will be conditioned on K, the number of 
simultaneously transmitting nodes. In general, K is a random 
variable, whose distribution is a function of the traffic char- 
acteristics, e.g, traffic load, traffic distribution, transmission 
control scheme, etc. In the simple case in which each node 
transmits with identical probability P t , K has a binomial 
distribution. Once the distribution of K is given then we can 
determine the SEP as P s = J2k=i p ( K ) p s( K )- 

From (O, the received signal at the destination q m is 

N 

y(fan;m) = ^ m 6 m ^ \a ml \ 2 s m 

i=l 

N K 

i=l 3=1 

where the first term is the desired signal and the remaining 
terms represent interference and noise. Recall that aji ~ 
CM (0, (Jq J . Since Sj is a PSK symbol, the magnitude of cijiSj 
is <j 2 s \a,ji\ and its phase is still uniformly distributed in [0, 2ir]. 
Thus, cijiSj ~ CM (0, u 2 a (jf). Therefore, 

K 

m = a ^ s J +m~CM (0, al) (11) 

3=1 
#™ 

where^A (K _i) a ^ + C7 2 

Given a m i, the instantaneous SINR, 7, equals 

a 2 b 2 la -\ 2 \ 2 a 2 u 2 b 2 £ 2 <r 2 

where £ = J2i=i \ a mi\ 2 - 



Note that fi rn is of order l/N. As N — > oo, fj^b^a 2 — > 0, 
and 7 reduces to /J^^V,/^, which corresponds to the 
scenario of additive white Gaussian noise (AWGN). Thus, 
under certain transmit powers, no matter how large N is, the 
SEP of the proposed scheme is always lower bounded by the 
SEP under AWGN. 

Since |a m ,| is Rayleigh distributed, £ ~ Erlang(iV, cr 2 ). 
The pdf of the Erlang distribution is 

Erlang(M): /fa M) = fffiljp X ^ " (13) 
The moment generating function (MGF) of 7 is 



17) in we obtain 

Ps(K) = - 



+ — — )-^^ 



sm ^3 



N—l 



M-l 1 / c r/ ?r , „ / 2n\ 1 

-\/t— { 9 + tan C E 

7T V 1 + C 2 i 



M it V 1 + c LV 2 ' " ^ V " / [4(l + c)l n 

71— v /J 



JV-l 



- 1 j=i { + Cj 



(19) 



A4 7 ( S ) 



exp(s7)/«(0 d C 



where 



and 



/ 7T 

— cot 

c \M 



exp( 



Sit 2 fo 2 t2 2 



(14) 



rp 



( a £-?>)4*[2(»-j) + l] 



(20) 



(21) 



based on which, the average SEP for M-PSK symbols is J5J 

(M-1W 



Ps{K) = - 



M 1 ( - 



sm' 



f -/O 



x- 



AT-1 



a^(iV-l)! 



exp 



sin 99 / 
sin 2 (7r/M) i^bU^l 



Sin 2 <yS Mm&jU ^ + ^ 



Recalling that M^(s) > M^(s) > 0, we have P S {K) > 
P S (K). The result of JT9l is an upper bound of the exact SEP 

of flSJ. 

An even simpler upper bound for P S (K) can be obtained 
based on Eq. (5A.76) of |0: 

/ \ ~ N 

P S {K) < n .,l_„„A (22) 



71/ 



sin 2 ( 



2^ (A/~l)7t^ 



(15) 



which for BPSK becomes 



* 2 c 



, -JV 



(23) 



Since there is no closed-form expression for A4 7 (s) or 
P S (K), in the following we will make some approximations 
to simplify the above expressions. 

A. A Simple Bound for SEP 

Let us fix an e > 0, and define £0 such that P(£ < £0) = e - 
Also, let us define 



2 ft 2 £2 2 



/i 2 6 2 (J 2 

C — C^C; 



Mm^CT 2 + CT 2 £/£ Mm^rnCT 2 +O- 2 /^0 



(16) 



When e is small, it holds with probability > 1 — e that 7 < 7. 
Since c 7 > and s is negative in the range of interest, we 
can always find a small enough e so that Ai^(s) > Ai 7 (s). 

Note that 7 ~ Erlang( N, a 2 a cx/) and thus the MGF of 7 is 
of the following simple form: 

My(a) = (1 -sa 2 a c^)~ N . (17) 

From dT3J, the SEP for M-PSK symbols based on 7 is 



1 Am-i)k/m sin 2 (^)cr 2 c^ 

P S (K) = - (1 ^ Sm UfKaC 7w 



(18) 



sm 99 



Defining c = sin (^)c 2 c^, and using the result of Eq. (5A. 



Remark: As /z m — ► 00, c 7 reduces to a 2 /a 2 , in which 
case the corresponding result of ( fT9l can be viewed as a 
lower bound when the transmit power of collaborating nodes 
approaches infinity. It shows that no matter how large the 
transmit power is, the SEP can never be smaller than this 
bound. The SEP floor is a result of the interference from other 
source nodes. To achieve lower SEP for a given K, one must 
increase N. Based on (1221 . this bound decreases approximately 
in a power-law fashion as N increases. 

V. Simulations 

In this section, we study the SEP performance of the 
proposed method via simulations, and also via the proposed 
analytical expressions. 

We assume the channels among nodes in a cluster are 
selected from zero-mean complex Gaussian processes, which 
are constant within one slot, but vary between slots. Let us 
define 7! = cr 2 cr 2 / 'a 2 ^, which represents the average SNR 
in the process of information sharing, and define 72 = 
N 2 /j, 2 n b 2 n a 2 a^/a 2 to represent the asymptotic average SNR 
(when N — > 00) at the receiver. Note that 72 is independent of 
N since fi m is of the order of l/N. Eq. (fT2l can be rewritten 
by 

e/N 2 



7 



+ 7 2 



(24) 



where £ = ~ Erlang(7V, 1). Then, the SINR is deter- 

mined only by 71, 72, K and N. Each packet contains BPSK 
symbols, so SEP is equivalent to BER. We take e = 0.01. 
Also, we assume perfect knowledge of channels, number of 
source nodes and destination information. Only one beampat- 
tern is formed in each slot. For simulation-based BER, we 
perform a Monte-Carlo experiment consisting of 10 6 repeated 
independent trials. 

Fig. Q] shows the BER versus 72 estimated from the network 
simulation (o line) when K = 4 nodes transmit all the time. 
The parameter 71 is fixed at 20 dB. The estimated BER is in 
perfect agreement with the analytical result for the exact SEP 
of (TTBTlCy line); in fact the two lines are indistinguishable. 
The upper bound on the exact SEP, computed by ([T9T l, is shown 
as the solid line. One can see that e = 0.01 can guarantee a 
tight bound under various parameters and SNR ranges. The 
simple upper bound computed via (f23l > is also shown (dashed 
lines). 

Extensive simulations confirm that the simulation-based 
BER and analytical SEP match well under a wide variety 
of scenarios. Thus, in the following we will simply use the 
analytical result of (TT~5T > to study the performance of the 
proposed method. 

Fig. |2] shows how the BER depends on the number of 
collaborating nodes for 71 = 20 dB and different values of 
72. Fig. [3] shows how K affects BER, where 71 = 72 = 20 
dB. As K increases, the SEP increases. Fig.|4]shows how BER 
changes with 71, where 72 = 20 dB and K — 4. Recall that 
cri = cr^a^(K — 1 + 7-f 1 ). K plays a dominant role in the 
interference (when K > 1). As observed in Fig. |4] the SEP 
decreases only slightly with the increase of 71. 

VI. Conclusions 

We have proposed a scheme for wireless ad hoc networks 
that uses the idea of collaborative beamforming and at the 
same time reduces the time needed for information sharing 
during the collaborative phase. We have provided an analysis 
of the SEP, which shows how the performance depends on the 
number of collaborating nodes, the number of simultaneously 
source users and noise levels at collaborating nodes and the 
final destination node. 
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Fig. 1. BER vs. 72 (K = 4, 71 = 20 dB); N = 8, 16, 32; empirical results, 
analytical exact results and upper bounds. 
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Fig. 2. BER vs. N (K = 4, 71 = 20 dB); N = 8, 16, 32; analytical exact 
results and upper bounds. 




exact, analytical 

the upper bound, e=0.01 

the bound of (23) 



Fig. 3. BER vs. K (71 = 72 = 20 dB); N = 8, 16, 32; analytical exact 
results and upper bounds. 
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Fig. 4. BER vs. 71 (K = 4, 72 = 20dB); TV = 8, 16, 32; analytical exact 
results and upper bounds. 



